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ABSTRACT 

We construct a filame nt catalogue using an extension of the halo based filament finder 
of lZhang etall(l200l . in a 250 Mpc/h side N-body simulation, and study the proper- 
ties of filaments ending upon or surrounding galaxy clusters (within 10 Mpc/h). In this 
region, the majority of filamentary mass, halo mass, and galaxy richness centered upon 
the cluster tends to lie in sheets, which are not always coincident. Fixing a sheet width 
of 3 Mpc/h for definiteness, we find the sheet orientations and (connected) filamentary 
mass, halo mass and richness fractions relative to the surrounding sphere. Filaments 
usually have one or more endpoints outside the sheet determined by filament or halo 
mass or richness, with at least one having a large probability to be aligned with the 
perpendicular of the plane. Scatter in mock cluster mass measurements, for several 
observables, is often correlated with the observational direction relative to these local 
sheets, most often for richness and weak lensing, somewhat less for Compton decre- 
ment, and least often for velocity dispersions. The long axis of the cluster also tends 
to lie in the sheet and its orientiation relative to line of sight also correlates with mass 
scatter. 
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for our purposes, we obtain a filament catalogue, and 
consider those filaments connected to or in the vicin- 
ity of galaxy clu sters. Our work is most clo s ely re - 
lated to that of IColberg. Krughoff fc Connolly! l|2005h : 
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estimates of th eir multiplicity (|Lee1l2006l ; IShen et al.ll2006l) . 
anisotropy (e.g lLee fc Springell (120091)) ? the merger rates 
of halos into them ( Song fc Led 20100 and properties in 
non-G aussian 
l2010h . 

Galaxy clusters (dark matter halos with mass 
M > lO 14 /^ 1 Mq) are of great interest for many reasons, in 
part because of their sensitivity to cosmological parameters, 
but also as hosts of the most massive galaxies in the uni- 
verse, as environments for galaxy evolution more generally, 
and as the largest virialized objects in the universe with cor- 
respondingly special astro - physic al processes and histories 
(for a review see e.g. IVoitJ (|2005l )). Galaxy clusters tend to 
lie at nodes of the cosmic web, with matter streaming into 
them fr om filaments fe.g. Ivan Haarlem fc van de Wevgaert! 
i| 19931 ); iDiaferio fc Gellerl (|l997l ); IColberg et all (|l99gh yT 
Although the universe is isotropic and homogeneous on 
large scales, around any individual cluster there will be 
directionally dependent density fluctuations due to the 
condensation of filamentary and sheetlike matter around it. 
Our interest here is in characterizing this nearby (within 
10 Mpc/h) filamentary environment of galaxy clusters. 
This environment feeds galaxy clusters and is also un- 
avoidably included for many observations of the cluster 
at its center. This correlated environment is one source of 
the observationally well known "projection effects," which 
have plagued optical cluster finding star ti ng with Abel] 
1958|) and later (e.g. iDalton et all (119921); iLumsden et a' 



I Aragon-Calvo. van de Wevgaert fc Joned (|2010h - They used 
simulations to measure counts of filaments (found via dif- 
ferent algorithms) ending upon clusters and average fila- 
mentary profiles and curvature. We go beyond these to 
measure the statistics of the local geometry of filaments 
around their cluster endpoints. Related studies of fila- 
ment ge ometry, particularly for superclusters ar e foun d 
in e.g.. I Aragon-Calvo. van de Wevgaert fc Jonesl (|2010l) ; 
I Aragon-Calvo. Shandarin fc Szalavl (|201ol ). the former also 
discuss the tendency of filaments around voids and clusters 
to lie in sheets. We find that most of the filamentary (and 
halo) material in a 10 Mpc/h sphere around clusters lies 
in a plane, presumably the one from which the filaments 
collapsed, and investigate different ways of defining such a 
plane's orientation. Many measures of cluster masses include 
the cluster environment and as a result scatter the mass 
from its true value. In mock observations on simulations, we 
find that line of sight dependent scatter in measured cluster 
masses, for several methods, is often correlated with the an- 
gle between the line of sight and these locally defined planes. 

In §2 we describe the simulations, mock observations, 
and filament finder. In §3 we describe the statistical proper- 
ties of the filaments and matter distribution around clusters, 
in §4 we consider the geometry of the filament, mass and 
richness distributions within 10 Mpc/h of each cluster, fo- 
cussing particularly on planes maximizing these quantities, 
in §5 we compare scatter in cluster masses to orientation of 
observations with these planes, and in §6 we conclude. 



19921 ); Ivan Haarlem et all (|l997l); White et all Jl999l)) 



weak leasing, e.g., (Reblinskv fc Bartelmannl 2 SIMULATIONS AND METHODS 



clust er - 
19991; iMetzler. White fc Lokenl BoOll; iHoekstral 



de Putter fc White 1200 
Becker & Kra vtsov |2 01 
l|Sunvaev fc Zel'dovich 



iMeneghetti et al 



2001 



2010 



cl uster 
19721. Il980l) 



Sunyaev-Zel'dovich 
(SZ) flux mea- 



surements, 



"dWhite. Hernquist fc Springell 



Holder. McCarthy fc Babull 120071 ; lHallman et all 



:.g.. 



2002 



2007 



Shaw. Holder fc Bodd l200sl) and cluste r velocity disper- 



dCenl 1 19971; iTormenl 1 19971 ; iKasun fc Evrardl 



120051 ; H iviano et all l200q )T The environments of clusters 
have been studied within several contexts and using 
several methods, e.g. galaxy and dark matter densit y 
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Here we describe our findings on local cluster en- 
vironments obtai ned by implementin g the halo-based fil- 
ament finder of IZhang et all ((2009) in a high resolu- 
tion N-body simulation. After refining the finder slightly 



2.1 Simulation 

We use a dark matter only simulation, in a periodic 
box of side 250 M pc/h with 2 048 3 particles evolved us- 
ing the TREEPM (|Whitdl2002T ) code, and provided to us 
by Martin White. It is th e same simulation as used in 
IWhite, Cohn fc SmitJ l|2010l ) (hereafter WCS), which can be 
consulted for details beyond those found below. The back- 
ground cosmological parameters are h = 0.7, n = 0.95, 
£l m = 0.274, and <rg = 0.8, in accord with a large number 
of cosmological observations. The simulation has outputs 
at 45 times equally spaced in ln(a) from z = 10 to 0. We 
focussed on z = 0.1, in part to allow comparison with ob- 
servational quantities in §5. H alos are found usi ng a Friends 
of Friends (FoF) halo finder (|Davis et all 1 19851 ). with link- 
ing length b — 0.168 times the mean interparticle spacing. 
Masses quoted below are FoF masses. 

Resolved subhalos in this high resolution simula- 
tion are of importance for the observational compar- 
isons in §5, and for measurements of galaxy properties 
around and in the clusters. Subha los are found via FoF6d 
| Diemand. Kuhlen fc Ma dau 2006]) , with the specific imple- 
mentation as described in the appendix of WCS. The sub- 
halos correspond to galaxies with luminosities > 0.2L* at 
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z = O.lQ, and match observations as described in WCS. The 
halo and subhalo catalogues and dark matter particles can 
be combined to produce mock observations for six cluster 
mass measures. These are (see WCS for specifics and tests 
of the catalogue): tw o richnesses (one using the MaxBCG 
i|Koester et al.l l2007f > algorith m based upon colorfl and 
one b ased upon spectroscopy |Yang. Mo fc van den Boschl 
2008)), SZ flux or Compton decrement (flux within an an- 
nulus of radius rigof>, the radius within which the average 
mass is greater than or equal to 180 times background den- 
sity), weak lensing (using an SIS or NFW model to assume a 
cluster lens profile and then fitting for a velocity dispersion 
and then mass), and two velocity dispersions (one based on 
a simple 3 — a clipping, the other on a more complex method 
using phase space information to reject outliers and calcu- 
lating mass usi ng a measured harmoni c radiu s as well, based 
on me t hods of [d en Hartog & Katgcrt ( 1996); iBiviano et al.l 
<|2006h ; IWoitak et all (|2007l )); more detail is in WCS. We 
will use the mass measurements by WCS via these meth- 
ods, taking cylinders of radius rigob when a radius choice is 
required. Just as in that work, lines of sights for clusters are 
removed a more massive cluster has its center within this 
radius along the observational line of sight. 

2.2 Filament finder 

We find fi laments usin g an e xtension of the method de- 
scribed in Zh ang et al.l (|2009h . They identify filaments as 
bridges in dark matter halos above a threshold halo mass 
overdensity, of length up to 10 Mpc/h. It is analogous to 
the spherical overdensity finder for clusters, where the clus- 
ter radius is taken to be that where the average density 
around the central point drops below some threshold; here 
the filament radius is where the average density along the 
cylinder axis drops below some threshold. Just as there are 
many different halo finders, there is no unique filament finder 
or definition. This finder is but one of many different ones 
present in the literature, which are not only based upon such 
bridge-like definitions, but also include finders constructed 
around filtering procedures, potential or density gradients , 
dynamical information, and more (see IZhang et all l|2009l ) 
for some comparisons between their finder and others). Even 
for a given filament finder, catalogues often must be speci- 
fied by the finder parameters as well (e.g. smoothing length 
for density or potential based finders, unbinding criteria for 
dynami cally based finders, etc.). We use the parameters as 
given in I Zhang et~aH (120091). 

The algorithm of IZhang et"aH (|2009h is as follows: ha- 
los are ordered most to least massive. All halos with mass 
> 3 x 10 w h- 1 M Q are includeojf], mass in the following only 
refers to that in halos with this mass or above. Starting with 
the most massive halo ("node"), all halos within 10 Mpc/h 
but at least 3 Mpc/h away in radius (or r2ooc if greater) are 
considered as potential endpoints. For each potential end- 
point, the cylinder radius is varied, up to 3 Mpc/h, to get the 



highest overdensity of halo matter in the cylinder between 
the node and potential endpoint. This maximum density is 
then compared to a minimum overdensity (5 times back- 
ground matter density in halos), and if over this minimum, 
this endpoint and its radius are kept. If no potential end- 
points have a halo mass density for their filament greater 
than the minimum overdensity then the algorithm moves to 
the next node. Once all such maximal filaments are found 
for a given node, the filament with the largest density is 
kept. The filament is then truncated: its new endpoint is 
the most massive halo within it which has at least 3 other 
halos between in and the central node, and which is at least 
3 Mpc/h away from the central node. All filament members 
are then removed from the list of potential future filament 
members or endpoints around any node. The endpoints are 
not removed from the list of possible endpoints for other 
nodes, but are removed from the list of possible endpoints 
associated with this node. This procedure is repeated until 
no more new filaments are found around the node. 

As this procedure frequently produces many more fila- 
ments than were evident by eye around clusters (sometimes 
over thirty around a single cluster), we incorporate a growing 
and merging procedure as well. After finding the filaments 
of maximum density around a given node, we grow out the 
filament radii until the average mass density in halos within 
the cylinder stretching to the filament endpoint drops to less 
than the minimum overdensity, or the maximum 3 Mpc/h 
radius is reached. Halos lying in two or more such extended 
filaments are assigned to the one whose axis is closest. Fila- 
ment endpoints with length I and a perpendicular distance 
d± to another (longer) filament's axis such that d±/£ < 3/10 
(the maximum width/maximum length in the algorithm) are 
merged into the longer filament, unless the shorter filament's 
endpoint has other filaments extending out of it. (This al- 
lows filament radii >3 Mpc/h.) These new filaments are then 
given a central axis determined by the center of mass of the 
filament; filaments whose endpoints do not have additional 
filaments extending out of them and whose endpoints are 
within 25 degrees of each other are merged. This is done in 
order of closest to most distant pairs; if > 2 filaments are 
within this range, the two closest are merged, then centers 
of mass are recalculated to see if the remaining filaments are 
within the minimum distance, and so on. 

The resulting filaments are regions connecting halos 
with halo mass overdensity at least 5 times the background 
halo mass density, and which are less than 10 Mpc/h long. 
The full catalogue at z = 0.1 has ~ 30,000 filaments and 
~ 44, 000 endpoints, with 45% of the halo mass fraction in 
filaments and 36% of the halos (in number fraction) in fil- 
aments. 60% of the ~ 1.2 x 10 6 halos above the minimum 
mass cut are not either endpoints or in filaments, with the 
most massive of these having M — 2.6 x 1O 12 /i _1 M0.Q 

Several of the other finders produce filaments which 
can extend well beyond our 10 Mpc/h cutoff (e.g. 
IColberg. Krughoff fc Connolly! l|2005l ) found filaments out 



Approximately -18.5 in r band, see WCS for more discussion. 

2 Color assignments are estimated us ing the prescription of 
ISkibba fc Shethl J2009I) with evo l ution o f IConrov, Gunn fc White! 
ll2009l l: IConrov. White, fc GumJ feOlcMConrov fc Gumji 20101) . 

3 This is t he min i mum mass used by IZhang et alj J200£ ) con- 
verted (see I White! feOOlf) ') to our FoF definition. 



4 Analytic estimates of filamentary mass fractions mentioned 
above (which use other filament definitions) are not directly com- 
parable because the latter are based upon total mass; mass in 
halos above our minimum is only 40% of the mass in the box at 
z = 0.1. 
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to 50 Mpc/h, even longe r ones have been found by e.g, 
iGonzalez fc Padillal J2009I)), some have restrictions on fil- 
ament nodes (e.g. IColberg. Krughoff fc Connollvl (|2005h fil- 
aments end only on clusters). Our catalogue has straight 
filament segments < 10 Mpc/h in length, built out of dark 
matter halos above some minimum overdensity which em- 
anate from clusters and other endpoints. Longer filaments 
could presumably be constructed as chains of our shorter 
ones, augmented by a condition on how much a filament 
can bend before it is considered instead to be two sepa- 
rate filaments meeting at a node. The length restriction of 
our finder also affects breakdowns into mass fraction in fil- 
aments, nodes, and so on, as some of our nodes will instead 
be filament members if the filaments are extended this way. 

The work most similar to ours in focus, 
studying clusters as fila men t en dpoints, is in 
IColberg. Krughoff fc Connollvl (|2005l ). some re- 
lated results are also found and compa red in 
lAragon-Calvo. van de Weygaert fc Jones! (|201dh (see 



also Sousbie. Pichon fc Kawaharal 1 2010f ). who found a clus- 
ter as an intersection of filam e nts in observational data). 
IColberg. Krughoff fc Connollvl (|2005T l found filaments by 
looking for matter overdensities by eye between cluster 
endpoints and measured a wide range of filament statistics, 
including the number of filaments per cluster as a function 
of mass, stacked filament profiles, length distributions, 
and the fractions of cluster pairs conne cte d by filaments. 
lAragon-Calvo. van de Weygaert fc Jones! ([2010) found 
filaments using a Multiscale Morphology Filter (see their 
paper for details) and consi de red s imilar quantities to 
IColberg. Krughoff fc Connollvl (2005), and in addition 
introduced a classification for filaments. 



3 STATISTICS OF FILAMENTS AROUND 
CLUSTERS 

Our finder is well suited to characterize the local environ- 
ment of clusters, our target of study here. Of the 242 clus- 
ters (M > 1O 14 /i _1 M0) in our box, 226 are also nodes, with 
~ 1700 filaments. We restrict to these clusters below. The 
other 7% (16) of the clusters are within filaments linking 
two more massive clusters, in addition, 29 of the clusters 
have a cluster within a filament, and 41 cluster pairs are 
within a 10 Mpc/h radius of each other. We use the term 
"connected" filamentary mass to refer to halo mass within 
a filament connected directly to a cluster, up to and includ- 
ing its other endpoint.0 In addition to connected filaments 
around a cluster, within the 10 Mpc/h sphere we will also 



5 The finder, even with modifications, still produced some con- 
figurations which we modified with post-processing. For example, 
sometimes a filament would be found with a large "gap" in the 
center, where the gap is due to a previously found filament be- 
tween two other clusters which crosses the region. Even with this 
gap, the new filament is above our overdensity threshold. As the 
previous and new filaments seem joined and perhaps one object, 
we added all the mass (within 10 Mpc/h) of any previously found 
filament which came within 3Mpc/h to the connected filamen- 
tary mass of the cluster, this happened for < 10 of our ~7000 
filaments. 
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Figure 1. Top: distribution of number of filaments per cluster 
(halos with M > 10 14 h~ 'Mq). Bottom: number of filaments as 
a function of mass for all halos which are filament endpoints. 



consider all filaments and their endpoints, all halos above 
our minimum mass of 3 x 10 10 h^ 1 Mq, and all galaxies. 

In the 10 Mpc/h spheres surrounding clusters, con- 
nected filaments constitute ~ 70% of the halo mass on av- 
erage, but with a very broad distribution of values for indi- 
vidual clusters. A line passing through the 10 Mpc/h shell 
centered on a cluster will hit one of the original connected fil- 
ament cores (from the first step of our algorithm) about 10% 
of the time on average, and one of the grown and merged 
filaments closer to ~ 30% of the time, with a wide spread 
as well. All (not only connected) filaments in this sphere 
contain closer to ~ 90% of the halo mass, with much less 
cluster to cluster scatter. (These additional filaments con- 
nect two other halos, which may or may not lie within the 
sphere themselves, rather the cluster and another halo.) In 
10 Mpc/h spheres around 10,000 random points, in compar- 
ison, the filaments have a halo mass fraction ranging from 
60% to 95%. 

The distribution of number of connected fila- 
ments around clusters, with our finder, is shown at 
the top of Fig. [T] clusters tend to have 7-9 fila- 
ments. We find more massive halos have more fila- 
ments ending upo n them, shown in Fig. [T] b ottom , 
just as found bv IColberg. Krughoff fc Connollvl l|2005h : 
lAragon-Calvo. van de Weygaert fc Jones! (|20101 i with their 
different finders. In addition, connected filaments around 
clusters tend to be shorter than their counterparts for the 
much less massive nodes. 

The large number of filaments found by the algorithm 
can be compared to a simplified picture where nodes are fed 
by a small number of filaments (e.g. 3 or fewer, iKeres et al.l 
l|2005l . 120091 ): iDekel et ail (|2009h ). The mass fraction in the 
largest 2 or 3 filaments is substantial, leading to a partial 
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corresponding to the largest eigenvalue) tends to align with 
the correspo nding direction of the fi l amentary inertia ten - 
sor. (See also Ivan de WevgaertJ (l2006l);lHahn et all (l2007tj); 
lAragon-Calvo et alj (l2007bl); |Paz. Stasvszvn fc Padillal 
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0.4 0.6 
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Figure 2. Cumulative fraction of filamentary mass in 1 (solid 
line), 2 (dashed line) and 3 (dotted line) most massive cluster fil- 
aments, as fraction of the number of clusters. For example, about 
half of the clusters have at least 60% of their connected filament 
mass in their two largest filaments and at least ~ 75% of their 
mass in their three largest filaments. 



reconciliation of these pictures, as seen in Fig. [2l that is 
about half of the clusters have at least ~ 75% of their mass 
in their three largest filaments. 

More massive halos have more filaments around them, 
more matter in filaments, and more matter around them 
generally, and although the number of filaments for clusters 
can be quite large, a significant fraction of the filamentary 
mass is found within the three largest filaments. 



4 PLANAR GEOMETRY AROUND CLUSTERS 

Filaments provide an anisotropic environment for galaxy 
clusters. Some approximate trends in the filamentary distri- 
bution are accessible via the inertia tensor of its mass, even 
though filaments are not expected to fill out an ellipsoid. 
For our clusters, the moment of inertia tensors tend to 
have two relatively large eigenvalues and one smaller one 
(corresponding to axis ratios a > b ~ c, the classic prolate 
cluster shape, there are ma ny studies of cluster ellipticities, 
see e.g. IJing fc Sutd (|2002l U In comparison, for connected 
filaments attached to our clusters, the middle eigenvalue of 
the inertia tensor tends to be smaller so that the filament 
distribution is "flatter" than the cluster it surroundsQ 
The long axis of the cluster has a tendency to lie within 
the "flat" directions of the filamentary distribution, and 
the eigenvector of the cluster's inertia tensor that is per- 
pendicular to the long and middle axes of the cluster (i.e. 



6 The connected filament distribution becomes more and more 
cylindrical with decreasing (well below 10 14 h~ 1 Mq) central halo 
mass, with the two largest eigenvalues tending to become equal, 
and the third becoming smaller and smaller. One reason is that 
lower mass halos are expected to be within filaments, rather than 
to serve as endpoints; the algorithm used here will tend to break 
these longer filaments up into more segments as mentioned earlier. 



1 20081 ): lAragon-Calvo. van de Weygaert fc Jones! (|201(t ); 
as our nodes will sometimes be members of filaments in 
other finders, some of these alignments are relevant filament 
member alignment discussed therein.) 

Visual inspection of many of our clusters sug- 
gests that the majority of their filamentary mass 
lies within sheetline regions, presumably those from 
which they condense d (see for example some cases 
illustr ated in lAragon-Calvo. Shandarin fc Szalavl 
(2010 ), and discussion of different fi l ament types in 
lAragon-Calvo. van de Wevgaert fc Jones! l|2010h and their 
"grid" and "star" configurations). 

To quantify this planarity, we consider four definitions 
of planes, regions extending ±1.5 Mpc/h above and below 
the central cluster and out to the edge of the local 10 Mpc/h 
sphere. We choose their orientations (normals) so the planes 
contain the maximum of either 1) connected filament mass, 
with extra constraints described below, 2) all filamentary 
mass including endpoints, 3) total halo mass, or 4) number 
of galaxies, within the 10 Mpc/h sphere. The connected fila- 
ment mass plane has its normal chosen to be perpendicular 
to the axes constructed out of a pair filament endpoints; this 
definition has stronger correlations with observables (dis- 
cussed later) than using pairs of connected filaments without 
their endpoints, or using the plane maximizing connected fil- 
ament mass with no other constraints. The mass in the plane 
(or richness, when using galaxies) does not include that of 
the central cluster, as our interest is in the cluster's envi- 
ronment. In Fig. [3]the objects used for these four choices of 
plane are shown for a cluster of mass 2.7 x 10 14 M Q //i. It has 
about 84% of its mass in the connected filament plane. 

These four planes tend to have similar orientations, with 
the all filament and halo mass planes are most often aligned 
(over 96% clusters have these two normals within 30 de- 
grees). This is not surprising given the dominance of fila- 
mentary mass in the 10 Mpc/h sphere around the cluster 
noted earlier. For a given cluster, the largest misalignment 
between any pairs of planes tends to be between its con- 
nected filament plane and one of the other planes, which for 
15% of the clusters differs by another plane by more than 60 
degrees. For most clusters it thus seems that the connected 
filaments are not as closely aligned with the other planes, 
which extend further out into the sphere. Plane pairs besides 
the closely aligned all filament and halo mass plane have on 
average 5 — 10% of the clusters mismatching by >60 degrees. 

The mass or richness fractions in these planes is sig- 
nificantly higher than the fraction (~l/5) of volume which 
the plane occupies in the sphere. The distribution of con- 
nected and total filament mass fractions, in the correspond- 
ing planes, for our clusters is shown at top in Fig. 21 at bot- 
tom is the distribution for the total halo mass plane. Also 
shown at bottom is the mass fraction for halo mass planes 
constructed around 10,000 random points (rescaled to have 
the same area under the curve), which is smaller on average 
than around the clusters. The richness fraction, not shown, 
peaks slightly more sharply than the halo mass fraction, but 
at a lower fraction (~ 60%). For all plane definitions, 80% 
of the clusters have more than 60% of their mass (or 55% of 



© 2010 RAS, MNRAS 000, 000-000 



6 Noh and Cohn 



connected filaments 



all filaments 




x 5 
total halo mass 





u 

CD 


40 


CO 




3 


30 








O 


20 


U 




CD 




-Q 


10 


G 




3 
G 






- filamentary mass 
: conn fil 


fraction n 


1 ' 1 ' 
1 1 






I all fil 






F- 1 — i — 1 — 








l 

. J 






1 

L 




1 1 ' — 


1 




i 







0.4 0.6 0.8 1 

fraction of filamentary mass w/in 1.5 Mpc/h of plane 




0.4 0.6 0.8 1 

fraction of total mass w/in 1.5 Mpc/h of plane 



Figure 3. Four types of objects used in constructing planes in a 
10 Mpc/h radius sphere centered on a 2.7 X 10 14 h~ 1 Mq cluster. 
Left to right, top to bottom are halos in connected filaments, 
halos in all filaments, all halos above 3 X 10 10 /i -1 Mq mass cut, 
and galaxies above 0.2L* cut (galaxies in cluster are not shown). 
Point size is proportional to halo mass or, for richness, halo infall 
mass (which determines luminosity, see WCS). About 84% of the 
cluster's connected filament mass is in the connected filament 
plane. 

their richness) in these planes; about a quarter of the clus- 
ters fell below this fraction for at least one plane definition. 
Our choice of plane height, ± 1.5 Mpc/h to give 3 Mpc/h 
in total, was motivated by the characteristic scale of clus- 
ter radii. We explored mass plane heights from 1-3 Mpc/h 
(total plane widths 2-6 Mpc/h), and found that the total 
halo mass fraction scaled as M p i ano /M S ph erc ~ height 1 ' 4 . It 
would be interesting to understand this scaling in terms of 
intrinsic filament profiles. 

The clusters with large plane misalignments (by > 60 
degrees) have low mass or richness fractions, or larger mass 
within 3 Mpc/h of the normal of the connected filament 
plane (but outside of it) almost twice as often as in the full 
sample (i.e. in ~2/3 of the clusters with mismatched planes). 
The misaligned plane clusters have only slightly more often 
a recenlQ merger or a larger intrinsic cluster flatness (as 
measured by its inertia tensor), they were equally likely to 
have other clusters within 10 Mpc/h as in the full sample. 

The connected filament plane's normal, similar to its 
counterpart for the connected filament's inertia tensor, tends 
to be aligned with its counterpart for the cluster's mass in- 
ertia tensor, and the cluster's long axis is likely to lie in the 
filament plane. The cluster galaxy positions, have an iner- 
tia tensor (setting mass to one) which appears uncorrelated 
with this plane, but restricting to more luminous (> 0.4L*, 
see WCS for detail) galaxies gives an inertia tensor whose 
"most fiat" (perpendicular to eigenvector for largest eigen- 
value) direction prefers alignment with the normal to the 
connected filamentary plane, and whose "long" axis tends to 

7 Specifically, a satellite which has fallen into the cluster within 
the last time step, ~ 600Af yrs, which had at the earlier time at 
least 1/10 of the cluster's final mass at z = 0.1. 



Figure 4. Top: Fraction of connected filament mass in connected 
filament plane (solid) and fraction of all filamentary mass within 
all filament plane (dashed), both in the fiducial 10 Mpc/h sphere 
around clusters. The normals to these planes are within 30 de- 
grees for ~80% of the clusters. Bottom: The fraction of total halo 
mass (above our 3 X X0 W Mq/H cutoff) in the mass plane around 
clusters (solid), and its counterpart around 10,000 random points 
(dashed, rescaled to have the same number volume as cluster his- 
tograms). A large fraction of the filamentary mass and total halo 
mass in 10 Mpc/h spheres around clusters resides within this pla- 
nar region containing ~20% of the volume. 

be within the filament plane. The galaxy velocity dispersions 
are larger along the connected filament pl ane than perpen- 
dicula r to it, but are more correlate d (e.g. iKasun fc Evrardl 
l|2005[ ); IWhite. Cohn fc SmitJ (|201Ch > with the inertia tensor 
of the cluster itself. 

Not all filaments lie in these planes. Filamentary mass 
can extend outside of the plane, as mentioned earlier, as can 
filament endpoints. The fraction of filament endpoints ly- 
ing outside the connected filament plane is shown in Fig. [S] 
note that this does not preclude a significant amount of the 
filament's mass lying within the plane. There is also an in- 
creased likelihood for at least one endpoint to lie perpendic- 
ular to the connected plane, as shown in Fig. [5] bottom. The 
distributions in Fig.[5]are similar for the other plane choices. 
About 1/10 of the clusters have more than ~3% of their 
connected filamentary mass within a 3 Mpc/h radius of the 
normal to their plane but above or below the plane itself, 
which we refer to as perpendicular filaments below. In ad- 
dition 10 clusters have over 15% of their mass in a region 
within 6 Mpc/h radius of the normal, but outside the con- 
nected plane. 

As noted earlier (Fig. \%$, the two most massive con- 
nected filaments often do possess a large fraction of the con- 
nected filament mass. The plane defined by these two fila- 
ments coincides with the connected filamentary mass plane 
almost half of the time[f|For 1/3 of the clusters, however, 
less than half of the connected filament planar mass comes 
from these two most massive segments. So although the two 



We thank G. Jungman for asking us to measure this. 
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Figure 5. Top: fraction of filament endpoints lying outside of 
connected filament plane- many filaments do not have their end- 
points in this plane, even though a large fraction of mass is in 
this plane (see Fig. |4f. Bottom: angle to normal of connected 
mass plane, for filament closest to the normal; at least one fila- 
ment tends to be perpendicular to this plane. The corresponding 
distributions for other planes are similar. 



most massive segments have a preponderance of filamentary 
mass (Fig. [2]), their large mass is not wholly responsible for 
the dominance of planar structure. 

The persistence of the locally defined planes to larger 
radii can be studied by fixing the plane height and orienta- 
tion, and extending the plane out into a region of 20 Mpc/h 
in radius, and calculating the fractional mass in this larger 
plane within the larger sphere. The plane volume fraction 
of the sphere volume drops by about one half compared to 
its value in the 10 Mpc/h sphere, but the (all) filamentary 
mass, halo mass and richness fractions in their respective 
planes drop by even more, by a factor of ~ 40%. There are 
filamentary, mass or richness planes in this larger sphere of 
the same ±1.5 Mpc/h width which have more of the fila- 
mentary, mass or richness in them (and usually more than 
1/2 the mass fraction of those defined within 10 Mpc/h). 
These 20 Mpc/h filament and mass planes differ from their 
counterparts at 10 Mpc/h by over 30 (60) degrees one half 
(one quarter) of the time, with slightly smaller fractions for 
the corresponding richness plane. 

We did not find a more useful measure of isotropy 
in the plane (i.e. in the angular direction), although 
the moment of inertia tensor can indicate how much 
the planar geometry tends to cylindrical (related ques- 
tion s have been explored when classifying fi laments, 
e.g. l|Arag6n-Calvo, van de Weygaert fc Jonesl l2010l ) note a 
"star" geometry for sets of filaments). One possible conse- 
quence of isotropy, or its lack, in the plane will be discussed 
in the next section on mass measurements. 

In summary, as has been known, the mass around clus- 
ters tends to lie in filaments, which themselves tend to lie 
within sheets. We have taken a set sheet width centered 
on the cluster and maximized different quantities (filament 
mass, connected filament mass, total halo mass and galaxy 



richness) within a 10 Mpc/h sphere around each cluster. 
The resulting planes are not always aligned: the all filament 
plane and all halo mass plane are most likely to be aligned, 
and the largest disagreement between planes for any cluster 
is most likely to be between the connected filament plane 
and another plane. The long axis of the cluster tends to lie 
in the plane as well. Often a perpendicular filament is also 
present relative to the plane, with others also partially ex- 
tending out of the sheet. The rough cartoon of the filament 
shape around clusters is a planar structure with a few fila- 
ments sticking out, with a tendency for at least one filament 
to be close to the plane's normal direction. 



5 CORRELATED MASS SCATTER WITH 
LOCAL FILAMENTARY PLANES 

There are observable consequences of the filaments sur- 
rounding galaxy clusters: most cluster observations, aside 
from X-ray 0, will tend to include some of the cluster envi- 
ronment as well as the cluster itself. We saw above that the 
majority of the clusters have a preferred direction in their 
local (10 Mpc/h radius) environments, with a large frac- 
tion of their surrounding (connected or all filamentary, or 
total halo) mass or richness lying in a 3 Mpc/h sheet. The 
relation of this local structure to observables can be stud- 
ied by using the mock observations described in WCS. In 
that work, cluster masses were measured along 96 lines of 
sight, using six methods mentioned earlier: two richnesses, 
Compton decrement, weak lensing, and two velocity disper- 
sions. For individual clusters, WCS found correlated out- 
liers in the mass-observable relation along different lines 
of sight. (It should be noted that Compton decrement and 
weak lensing both can have significant contamination from 
beyond the 250 Mpc/h path measured within the box, so 
correlation with the local environment is likely smaller than 
found in WCS and below.) Some connection with environ- 
ment or intrinsic properties is seen: for the 8% of cases 
where at least two observables had a large (> 50%) devi- 
ation in mass from that predicted by the mean relation, an 
excess of nearby galaxies from massive or less massive halos 
and/or substructure (as det ected by the Dressler-Shectman 
jDressler fe Shectmanll 19881 ) test) were found relative to the 
population without these outliers. 

The filamentary structures and mass planes, and the 
mass fraction in them, provide an additional characteriza- 
tion of individual cluster environments. The WCS mock ob- 
servations along the 96 lines of sight of each cluster can 
now be compared to |cos#|, where 9 is the angle between 
the line of sight and the normal to these planes. In addi- 
tion to the normal to four of the planes mentioned above 
(connected filament mass, filamentary mass, halo mass, and 
galaxy richness above 0.2L*), we also consider a fifth pre- 
ferred direction, the angle to the nearest filament, and in this 
case use |sin0fu| of this angle (i.e. | cos f9| of the associated 
normal to the nearest filament). 

The rough expectation is that a cluster's measured mass 
along the sheet with the most filamentary mass or total halo 

9 X-ray structure might have some correlation as well, inasmuch 
as X-ray substructure is related to filaments which provide the 
cluster's infalling material. 
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mass (| cos6>| ~ 0) will be larger than that along the plane's 
normal vector. This correlation is not expected to be per- 
fect, as there is often a filament close to the normal, and 
the fraction and distribution of mass in the plane can vary. 
In particular, the planes aren't necessarily completely filled, 
and some directions through this plane might not intersect 
large amounts of mass (i.e. there might be a lack of isotropy 
in the plane as mentioned earlier) . For planes which are not 
isotropically filled, one might thus expect a triangular distri- 
bution of mass prediction (on the x-axis) vs. | cos^l (on the 
y-axis): with low mass values for all | cos#|, and high mass 
values for small | cos#| (along the plane). In addition, planes 
were defined only within 10 Mpc/h of the cluster, or less, and 
interlopers are expected to have an effect sometimes at 10 
times or more of that distance. These factors suggest that 
the alignment of an observational direction with a sheet may 
not be highly noticeable in observations, even if most of the 
local (filamentary and/or halo) mass lies within this sheet. 

Even with these contraindications, for many clusters we 
found a strong correlation for many mass measures with the 
angle between the line of sight and the locally defined planes. 
These strong correlations are seen not only for both mea- 
sures of richness, which in principle are closely localized to 
the cluster, but also for weak lensing, and to a lesser extent, 
SZ. Correlations are both less frequent and less strong for ve- 
locity dispersions. We show an example of one cluster's mass 
scatter for the six observables in Fig. [U The measured mass 
is calculated using scaling from the mean mass-observable 
relation for clusters in the simulation with M > 10 14 /i _1 Mq, 
and its value is shown versus | cos 9\, where 9 is the angle be- 
tween the observational direction and the connected filament 
plane's normal. The six panels show two richnesses, SZ, weak 
lensing, and two velocity dispersions. This 2.7 x 10 14 /i _1 Mq 
cluster, with 9 filaments, exhibits strong correlations for all 
six measurements. It has 84% of its connected filament mass 
and 72% of its halo mass in the connected filament plane. 

Given the noisiness of the data, we are mostly interested 
in general qualitative trends for the full set of 226 cluster 
nodes. We estimate correlations for each cluster in two ways. 
One is to use the correlation coefficient for (log M, | cos 6\), or 
the truncated set of points by the procedure described below, 
if that gives a lower absolute value (i.e. weaker value) for the 
correlation coefficient. These are shown for our example in 
Fig. [6]above. By eye, a correlation of < —0.25 appears to be 
a strong correlation, between —0.25 and 0.25 is often (not 
always) extremely noisy, and a correlation > 0.25 indicates 
an (unexpectedly) positive correlation. We use this division 
hereon. A positive correlation is unexpected as this means 
that measured cluster mass increases as the line of sight 
intersects less of the preferred plane. 

The distributions of these correlation coefficients, for 
the respective mass measurements in Fig. [6] and the con- 
nected filament plane, are shown in Fig. for all 226 cluster 
nodes. Also printed are the number of clusters with strong 
(negative), noisy and positive correlations for each measure- 
ment. The results are similar for all 5 choices of plane within 
the considerable noise The fraction of clusters having 



10 Relative to the connected filament plane shown, all other 
planes have more strong negative correlations for the two rich- 
ness based masses; for planes besides the plane perpendicular to 
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Figure 6. Example of mass scatter correlations: each point is a 
mass measurement for the same cluster along one of ~ 96 lines of 
sight, having angle 9 with the normal of the connected filament 
plane. The vertical line gives the true mass. The mass measure- 
ments are based upon (left to right, top to bottom): red galaxy 
richness, richness based on phase space, Compton decrement, 
weak lensing, phase based velocity dispersions and 3-ct clipping 
velocity dispersions. The mass axes for each measurement vary to 
cover the range of masses found for that technique, note the scales 
differ. Envelopes are fit to truncated sets of these points, both us- 
ing a chi-squared fitting (dashed line) and a shortest perpendicu- 
lar distance to the envelope (solid line), as described in the text. 
Where the two severely disagree (e.g. lower left hand box), one 
or both fits are bad. The correlation coefficients between |cos0 
and log 10 M/h are shown at upper right, they are the smallest in 
absolute value of of those for all or truncated points., or the trun- 
cated points. The cluster has mass 2.7 X 10 14 M Q /h, 9 filaments, 
and about 84% of its connected filament mass in the connected 
filament plane. 



strong negative or positive correlations, split according to 
type of mass measurement, is shown in table [1] with ranges 
shown for the 5 choices of plane. The composite correlation 
of log A/true — logMprod for all the clusters with |cos(9| fol- 
lowed similar trends, with a strongest correlation coefficient 
for both richnesses, then weak lensing; velocity dispersions 
and SZ are all similarly low. The amount of correlation be- 
tween line of sight and normal to various planes is corre- 
lated to some extent with the fraction of mass or richness in 
these planes, as might be expected. There is also a correlation 
between the strength of correlations of (logAf, |cos#|) and 
the alignments between planes for each cluster (not surpris- 
ingly, this depends upon the pair of planes being considered 
and the plane used to defined 9). Considering multiwave- 



then nearest filament (which is lower), there are more strong cor- 
relations for weak lensing and Compton decrement, and similar 
numbers for velocity dispersions. The plane perpendicular to the 
nearest filament has fewer negative correlations for weak lensing 
and Compton decrement and many fewer for velocity dispersions. 
For all 6 mass measurements, the median correlation for the plane 
perpendicular to the nearest filament is weaker (i.e. more posi- 
tive) than for the other four planes, by more than the scatter 
between the median correlations for other four. 
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property 


red richness 


phase rich 


SZ 


weak lensing 


phase v 


3<r v 


corrln < -0.25 


70-80% 


85-90% 


35-50% 


55-75% 


20-40% 


25-40% 


corrln < -0.25 or 
neg rh slope 


75-85% 


85-95% 


40-55% 


65-75% 


25-40% 


30-40% 


inertia corrln <-0.25 


80-85% 


75-80% 


40-50% 


>90% 


35-45% 


55% 


corrln > 0.25 


<3% 


<3% 


<2% 


1-6% 


1-5% 


2-7% 


corrln > 0.25 or 
large rh pos slope 


1-4% 


2-5% 


< 3% 


4-10% 


4-9% 


2-7% 


ill defined slopes 


5% 


~ 


~ 


~ 


45-50% 


50-60% 



Table 1. Cluster fractions with strongly negative (expected, rounded to the nearest 5%) and positive (unexpected, not rounded) 
correlation coefficients between |cos0| and measured mass, by observable, for filament planes, and effect of also considering the inverse 
slope of the right hand ("rh") envelopes of the (logM, | cos£>|) relation, when either strongly negative (expected) or positive (unexpected), 
and for directions associated with cluster inertia tensor. The range of values encompass those for planes defined with connected filament 
mass, all filamentary mass within 10 Mpc/h sphere, all halo mass within 10 Mpc/h sphere, galaxy richness, and the plane whose normal is 
perpendicular to nearest filament to line of sight. Also shown for planes are fractions of clusters with badly defined envelopes ("ill-defined 
slopes"— suggesting no correlation). The range for strongly negative correlations with two directions of the inertia tensor of the cluster 
itself (long axis of cluster, using sinS, or direction of eigenvector with largest eigenvalue) is also shown, see below in text. 
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Figure 7. Distributions of correlations between measured mass 
and |cos0| for the 226 cluster nodes, for six observables. Here 
is the angle between the line of sight and the normal to the con- 
nected filament plane. The correlation for each cluster is taken to 
be the one which is minimum in absolute value for all points or the 
truncated (as described in the text) set of points. The mass mea- 
surement methods are as in Fig. [6] i.e., left, right, top to bottom 
are red galaxy richness, phase space richness, Compton decre- 
ment, weak lensing, velocity dispersion using spatial information 
and 3cr clipping velocity dispersion. Also printed for each method 
are (left) the number of clusters with correlation < —0.25, (mid- 
dle) the number of clusters where the correlation's absolute value 
is less than 0.25 (and thus possibly noise), and (right) the num- 
ber of clusters where the correlation is > 0.25, i.e. both positive 
and large, indicating a higher mass estimate as the line of sight 
becomes more perpendicular to the maximal plane. The dashed 
vertical lines separate these three regions. 



length measurements together for each cluster, ~40-50% of 
the clusters have a strong negative correlation (i.e. the ex- 
pected sign) for at least 3 observables. 

For the planes, the correlation coefficient sometimes 



is low, even with a visible trend of measured mass versus 
cos#|. One apparent cause is the expected triangular enve- 
lope for the points described above. To identify this pattern, 
we considered slopes of approximate envelopes of the distri- 
butions, shown in Fig. [6] Points are binned in 8 approxi- 
mately equally filled 11 ! Icosr^l bins, in each bin < 2 points 
are discarded at large or small log M if separated from their 
nearest neighbor by more than 6 times the median separa- 
tion in mass in that bin (or the minimum separation if the 
median is zero). This threw out many of the notable out- 
liers. It also sometimes threw out other points, in a binning 
dependent way, but the number these points is small and 
not a concern as we are interested in average overall prop- 
erties. Points within 3<j of the median log M are then kept 
within each | cosf? bin. Straight line envelopes were then fit 
to both ends of each bin, either by minimizing perpendicular 
distance to the envelope or minimizing the chi-squared (note 
of log M (| cos 6\)). Envelopes for both methods are shown in 
Fig. [5] the cases shown where they strongly disagree corre- 
spond to one or both envelopes having bad fits. From hereon 
we restrict to envelopes based upon minimizing the perpen- 
dicular distance to the envelope. The resulting right and left 
hand inverse slopes are correlated with the correlation co- 
efficents of (log M, | cos 9\). We explored adding clusters to 
the negative (or positive) slope sample which have inverse 
slopes less than (or more than) the mean value of inverse 
slope for our correlation coefficient cutoff ±0.25; the small 
effect can be seen in table[T] (Sometimes the mean value had 
the wrong sign, e.g. for velocity dispersions for some choice 
of plane, which have large scatter, in this case the cutoff 
was set to zero.) We strove to be conservative in claiming a 
correlation, so that our estimates for the strength of these 
planar orientational effects tend to be lower bounds. 

Most of the time the envelopes found by our algorithm 
are reasonable to the eye, but sometimes they fail catas- 
trophically, and were caught by the goodness of fit estimator. 



11 As mentioned earlier, lines of sight where a more massive clus- 
ter is present within rigob are discarded. 
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The catastrophic failures seem to occur when no correlation 
is apparent between | cos#| and measured mass, as does an 
envelope close to vertical (inverse slope close to zero). The 
goodness of fits are the worst for the velocity dispersions, 
which have close to half of the clusters not allowing good fits 
for either the left or right envelopes, even when the good- 
ness of fit passes threshold, the envelopes are often close to 
vertical: i.e. the minimum or maximum velocity dispersion 
mass is similar either perpendicular to the maximum plane 
or looking through it. 

For the unexpected positive correlations, a positive in- 
verse left envelope slope can be understood by looking down 
a filament near the perpendicular to the plane (small angle, 
large mass) and then catching a "gap" in the plane (large 
angle, small mass). It is more difficult to understand > 0.25 
correlations or large positive right envelope inverse slopes 
(i.e. the largest measured mass closer to the perpendicular 
to the dominant plane). These do not dominate but are not 
uncommon: for any choice of plane, ~10-20% of the clusters 
have at least one observable with strongly positive inverse 
right hand slope or correlation (almost half of these are due 
to velocity dispersions). This dropped to < 5% (down to 1% 
using the plane perpendicular to the nearest filament or rich- 
ness) when requiring clusters to have at least 3 observables 
with either right hand positive slope or correlation (most 
often weak lensing and both dispersion measurements). 

Restricting to correlations, which are a cleaner and 
more conservative measurement, there are 45 clusters with 
a positive correlation for at least one measurement (usually 
velocity dispersions). These clusters differ from the full sam- 
ple in having, twice as often as the latter, high fractions of 
perpendicular mass to the connected filament plane and/or 
some pair of planes misaligned by 60 degrees or a recent 
merger (as defined earlier). They also slightly more often 
have another massive cluster within 10 Mpc/h, low mass or 
richness fraction in some plane, or are more fiat (as mea- 
sured by its inertia tensor, smallest axis/middle axis < 0.6). 
Fewer than a quarter of the clusters with a positive corre- 
lation for at least one measurement don't have one of these 
factors present, and some of these are close to our cutoffs, 
e.g. have more mass within 6 Mpc/h (rather than 3 Mpc/h) 
to the perpendicular to the connected filament plane than 
most clusters, or planes mismatching by almost 60 degrees. 
The "unexplained" strongly positive correlations occur for 
weak lensing and velocity dispersion mass measurements!* 2 ! 
Given the complexity of the cosmic web, and the small re- 
gion we use to characterize the cluster's environment, it is 
to be expected that our simple cartoon description will not 
always correlate precisely with observables. 

Similar correlations can be calculated using two axes de- 
fined from the inertia tensor for the cluster itself: the "long" 
axis of the cluster, and the normal to the direction of the 
eigenvector for the largest eigenvalue of its inertia tensor 

12 The fewest cases of strongly positive inverse slope or corre- 
lation occur for the plane defined using the perpendicular to the 
nearest filament to line of sight, suggesting that filaments close to 
the line of sight might be the cause of positive correlations, but 
again positive correlations did not always occur for these configu- 
rations. However, the plane perpendicular to the nearest filament 
also gives the fewest (except for richness) strongly negative cor- 
relations, i.e. its correlations are weaker in general. 



(pointing orthogonal to the longest and middle axes of the 
cluster). As mentioned earlier, these directions are corre- 
lated with the planes, with the "long" cluster axis tending 
to lie within them and the latter direction tending to align 
with the plane normals. Compared to the five planes above, 
the median correlation with mass scatter is stronger for red 
galaxy richness, weak lensing and the two velocity disper- 
sions, is similar for SZ, and brackets that for phase richness 
(the long axis of the cluster always has the stronger corre- 
lation of the two). The strength of effect for the "long" axis 
of the cluster is likely due to not only a plane being com- 
pared to the line of sight, but a specific high density axis 
within that plane; almost 90% of the clusters have a strong 
negative correlation for at least 3 of the 6 observables. The 
fractions of strong negative correlations for these two direc- 
tions determined by the cluster inertia tensor is also shown 
in table [T] 

In summary, the mass scatter for richness, Compton 
decrement, weak lensing and velocity dispersion measures is 
often correlated with the angle to these planes (most for rich- 
nesses, and least for velocity dispersions). The correlations 
aren't perfect and can sometimes be weak, or even of the op- 
posite sign than expected. In the latter case it is often also 
true that the different dominant planes (mass, connected or 
all filamentary halos, and richness) are not well aligned, or 
that a large filament extends perpendicular to the connected 
filament plane. Besides being correlated with each other, the 
planes and the mass scatters are also correlated with axes 
of the cluster's inertia tensor. 



6 CONCLUSION 

After implementing a filament finder on an N-body simula- 
tion, we studied the resulting filamentary environment for 
the 226 nodes which are also clusters (M > W 14 M Q /h). Fil- 
aments tend to lie in sheets, presumably those from which 
they condensed, providing a highly anisotropic environment 
for the cluster at their center. Within a 10 Mpc/h sphere, we 
identified sheets of width 3 Mpc/h, centered on each cluster, 
which maximize either total mass, connected filament mass, 
all filamentary mass, or galaxy richness. The all filament and 
halo mass planes are most often closely aligned, while the 
connected filament plane tends to be within the pair of least 
aligned planes for the majority of clusters. The direction 
of the filamentary and mass planes persist slightly as the 10 
Mpc/h spheres are extended to 20 Mpc/h. We measured the 
correlation of mass measurement scatters with the direction 
of observation relative to these planes for mock observations 
of richness, Compton decrement, weak lensing and velocity 
dispersions, via correlation coefficients and fits to the en- 
velopes of the measurements. Often there is a strong cor- 
relation between measured mass and direction to the local 
plane, in spite of the relatively small region (10 Mpc/h ra- 
dius) used to define the plane (again, this correlation might 
be overestimated for Compton decrement and weak lensing, 
which both can have strong scatter from distances larger 
than our box size). Strong correlations are least likely for 
velocity dispersions, and fitting envelopes to their distribu- 
tion of | cos 9\ versus log M tend to fail badly. This is perhaps 
not surprising because our finder doesn't include dynamical 
information. Alignments of observational direction with two 
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of the axes of the inertia tensor of the cluster also results in 
strong correlations with measured mass scatter. 

How these planes and correlations with scatter extends 
to higher redshift depends upon how the finder extends to 
higher redsh i ft. Th is is a subtle question as the finder of 
I Zhang et al.l (|2009h has a built-in scale: a cutoff for mini- 
mum halo mass. A full analysis of appropriate generaliza- 
tions is beyond the scope of this paper; two natural possi- 
bilities, however, are to leave the minimum mass alone, or 
to choose a minimum mass so that the ratio of the number 
of halos to the number of clusters (107 at z = 0.5, 25 at 
z — 1.0) remains the same, which gives a minimum mass 
of 8.2 x lO lo /i _1 M for z = 0.5 and 3.0 x lO 11 ft" 1 M for 
z = 1.00 Choosing the latter case (and luminosity cut at 
0.2L*), most of the trends persist to these higher redshifts, 
although the total number of filaments in the box decreases. 
For 2 = 0.5 and z — 1.0, the planar mass fractions around 
clusters are close to unchanged. For all 3 redshifts, there 
is a slight drop in richness fraction in the richness plane 
as redshift increases, and the halo mass fraction in planes 
around random points appears to grow, so that by z = 1.0 
it is comparable to that around the 25 clusters in the box 
at z = 1.0. For correlations of plane directions with cluster 
observations, the statistics are very noisy for z = 1.0, but for 
z = 0.5, the fractions of clusters with strong (expected) neg- 
ative correlations of angle with plane and mass scatteil 14 !, as 
in table[TJ tend to either remain the same in range or slightly 
increase (velocity dispersions do decrease in one case), the 
number of clusters with at least three negative correlations 
are close to unchanged for three planes, dropping for the 
richness and nearest filament planes, and positive correla- 
tion fractions are about the same except for (an increase 
for) velocity dispersions. Large plane misalignments are less 
common, but clusters with misaligned planes still are more 
likely to have smaller mass fractions in the plane or more 
perpendicular mass than the full sample. 

It would be interesting to determine whether this gener- 
alization to higher redshift is appropriate and then to under- 
stand the results in terms of the evolution of the filamentary 
neighborhood of the clusters and the clusters within them. 

The correlations between mass scatter and angle of ob- 
servation with the planes (and inertia tensor of the cluster) 
rely upon three dimensional information available to us as 
simulators. It would be very interesting to find a way to 
make this source of mass bias more evident to observers, 
perhaps by using a filament finder based upon galaxies di- 
rectly (amongst those mentioned earlier), and seeing how 
well they trace these planes, or by combining multiwave- 
length measurements . In depth studies underway of cluster 
environments such as iLubin et al.l (|2009l ) would be excellent 
datasets to apply and refine such methods. 
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ter for Physics for providing the place and opportunity for 
us to meet and discuss. YN thanks the Santa Fe Cosmol- 
ogy School for the opportunity to present this work, and we 
both thank the participants there for many useful comments 
and questions, and S. Habib and K. Heitmann for support 
in order to be able to attend. Martin White's simulations, 
used in this paper, were performed at the National Energy 
Research Scientific Computing Center and the Laboratory 
Research Computing project at Lawrence Berkeley National 
Laboratory. 



REFERENCES 

Abell CO., 1958, ApJS, 3, 211 

Altay, G., Colberg, J. M., & Croft, R. A. C. 2006, MNRAS, 
370, 1422 

Aragon-Calvo, M. A., Jones, B. J. T., van de Weygaert, R., 
& van der Hulst, J. M., 2007a, A & A, 474, 315 

Aragon-Calvo, M. A., van de Weygaert, R., Jones, B. J. T., 
& van der Hulst, J. M., 2007b, ApJL, 655, L5 

Aragon-Calvo, M.A., Shandarin, S.F., Szalay, A., 2010, 
larXiv:1006.4178l 

Aragon-Calvo, M.A., van de Weygaert, R., Jones, B.J.T., 

2010, larXiv: 1007.07421 
Bailin, J., & Steinmetz, M., 2005, ApJ, 627, 647 
Barrow, J. D., Bhavsar, S. P., & Sonoda, D. H., 1985, MN- 
RAS, 216, 17 

Basilakos, S.; Plionis, M.; Yepes, G.; Gottlber, S.; Tur- 

chaninov, V., 2006, MNRAS, 365, 539 
Becker, M.R., Kravtsov, A.V., 2010. larXiv:1011. 16811 
Betancort-Rijo, J.E., Trujillo, I., 2009, arXiv:0912.1051 
Bharadway, S., Sahni, V., Sathyaprakash, B.S., Shadarin, 

S.F., 2000, ApJ 528, 21 
Biviano, A., Murante, G., Borgani, S., Diaferio, A., Dolag, 

K., Girardi, M., 2006, A& A, 456, 23 
Bond, N., Strauss, M., Cen, R., 2009. larXivr0 903.3601 
Bond, N.A., Strauss, M.A., Cen, R., 2010, arXiv: 1003T3237l 
Bond J. R., Kofman L., Pogosyan D., 1996, Nature, 380, 

603 

Cen R., 1997, ApJ, 485, 39 

Chambers, S.W., Melott, A.L., Miller, C.J., 2000, ApJ, 544, 
104 

Choi, E., Bond, N.A., Strauss, M.A., Coil, A.L., Davis, M., 

Willmer, C.N.A., 2010, MNRAS, 406, 320 
Colberg, J., White, S.D.M., Jenkins, A., Pearce, F.R., 1999, 

MNRAS, 308, 593 
Colberg, J., Krughoff, K.S., Connolly, A.J., 2005, MNRAS, 

359, 272 

Colberg, J. M., 2007, MNRAS, 375, 337 
Colombi, S., Pogosyan, D., Souradeep, T., 2000, PRL 85, 
5515 

Conroy C, Gunn J.E., White M., 2009, ApJ, 699, 486 
Conroy C, White M., Gunn J.E., 2010, ApJ, 708, 58 
Conroy C, Gunn J.E., 2010, ApJ, 712, 833 
Costa-Duarte, M.V., Sodre Jr., L., Durret, F., 2010, 
larXiv:1010.0981l 

Dalton, G.B., Efstathiou, G., Maddox, S.J., & Sutherland, 

W.J., 1992, ApJL 390, LI 
Davis M., Efstathiou G., Frenk C.S., White S.D.M., 1985, 

ApJ, 292, 371 
Dekel, A., et al, 2009, Nature, 457, 451 



© 2010 RAS, MNRAS 000, 000-000 



12 Noh and Cohn 



den Hartog, R., Katgert, P., 1996, MNRAS, 279, 349 
Diaferio, A., Geller, M.J., 1997, ApJ, 481, 633 
Diemand J., Kuhlen M., Madau P., 2006, ApJ, 649, 1 
Dolag, K., Meneghetti, M., Moscardini, L., Rasia, E., & 

Bonaldi, A., 2006, MNRAS, 370, 656 
Dressier A., Shectman S.A., 1988, AJ, 95, 985 
Einasto, J., Klypin, A. A., Saar, E., Shandarin, S.F., 1984, 

MNRAS, 206, 529 
Faltenbacher, A., Gottlober, S., Kerscher, M., Muller, V., 

2002, A & A, 395, 1 
Faltenbacher, A., Allgood, B., Gottlober, S., Yepes, G., & 

Hoffman, Y., 2005, MNRAS, 362, 1099 
Faltenbacher, A., Li, C, Mao, S., van den Bosch, F.C., 

Yang, X., Jing, Y.P., Pasquali, A., Mo, H.J., 2007, ApJL, 

662, 71 

Feix, M., Xu, D., Shan, H., Famaey, B., Limousin, M., 
Zhao, H., & Taylor, A., 2008, ApJ, 682, 711 

Forero- Romero J.E., Hoffman Y., Gottlober S., Klypin, A., 
Yepes G., 2009, MNRAS, 396, 1815 

Gal, R.R., Lemaux, B.C., Lubin, L.M., Kocevski, D., 
Squires, G.K., 2008, ApJ, 684, 933 

Gay, C, Pichon, C, Le Borgne, D., Teyssier, R., Sousbie, 
T., & Devriendt, J., 2009. larXiv:0910.1728l 

Genovese, C.R., Perone-Pacifico, M., Verdinelli, L, Wasser- 
man, L., 2010. larXiv:1003.5536l 

Gonzalez, R.E., Padilla, N., 2009, larXiv:0912.0006l MN- 
RAS to appear 

Hahn, O., Porciani, C, Carollo, C. M., Dekel, A., 2007a, 
MNRAS, 375, 489 

Hahn, O., Carollo, C. M., Poricani, C, Dekel, A., 2007b, 
MNRAS, 381, 4 

Hahn, O., Teyssier, R., Carollo, CM., 2010, 
larXiv: 1002.19641 MNRAS to appear 

Hallman, E.J., O'Shea, B.W., Burns, J.O., Norman, M.L., 

Harkness, R., Wagner, R., 2007, ApJ 671, 27 
Hoekstra H., 2001, A&A, 370, 743 

Holder, G.P., McCarthy, I.G., Babul, A., 2007, MNRAS 
382, 1697 

Hopkins, P.F., Bahcall, N., Bode, N., 2004, ApJ 618, 1 
Jing, Y.P., Suto, Y., 2002, ApJ, 574, 538 
Jones, B.J.T., van de Weygaert, R., Aragon-Calvo, M.A., 
larXiv:1001.4479l 

Kartaltepe, J. S., Ebeling, H., Ma, C. J., & Donovan, D., 

2008, MNRAS, 389, 1240 
Kasun S.F., Evrard A.E., 2005, ApJ, 629, 781 
Keres, D., Katz, Weinberg, D. H. & Dave, R., 2005, MN- 
RAS, 363, 2 

Keres, D., Katz, N., Fardal, M., Dave, R., & Weinberg, 
D. H., 2009, MNRAS, 395, 160 

Koester B.P., et al., 2007, ApJ, 660, 221 

Lee, J., 2004, ApJL, 614, LI 

Lee, J., 2006, |astro-ph/0605697| 

Lee, J., Evrard, A.E., 2007, ApJ 657, 30 

Lee, J., Springel, V., Pen, U.-L., Lemson, C, 2008, MNRAS 
389, 1266 

Lee, J., Springel, V., 2009. larXiv:0911.4755l JCAP to ap- 
pear 

Lubin, L.M., Gal, R.R., Lemaux, B.C., Kocevski, D.D., 

Squires, G.K., 2009, AJ, 137, 4867 
Lumsden, S.L., Nichol, R.C., Collins, C.A., Guzzo, L., 1992, 

MNRAS 258,1 

Mead, J. M. G., King, L. J., & McCarthy, I. G., 2010, 



MNRAS, 401, 2257 
Mecke, K. R., Buchert, T., & Wagner, H., 1994, A&A, 
288, 697 

Meneghetti, M., Fedeli, C, Pace, F., Gottloeber, S., Yepes, 

G., 2010. larXiv:1003.4544l 
Metzler C, White M., Loken C, 2001, ApJ, 547, 560 
Murphy, D.N.A., Eke, V.R., Frenk, C.S., 2010, 

arXiv: 1010:2202 
Novikov, D., Colombi, S., & Dore, O., 2006, MNRAS, 366, 

1201 

Onuora, L.I., Thomas, P.A., 2000, MNRAS, 319, 614 
Pandey, B., Somnath, B., 2006, MNRAS, 372, 827 
Paz, D. J., Stasyszyn, F., & Padilla, N. D., 2008, MNRAS, 
389, 1127 

Pereira, M. J., Bryan, G. L., & Gill, S. P. D., 2008, ApJ, 
672, 825 

Pimbblet, K. A., 2005a, MNRAS, 358, 256 

Pimbblet, K. A., 2005b, Publications of the Astronomical 

Society of Australia, 22, 136 
Pimbblet, K. A., Drinkwater, M. J., & Hawkrigg, M. C, 

2004, MNRAS, 354, L61 
Poggianti, B.M., De Lucia, G., Varela, J., Aragon- 

Salamanca, A., Finn, R., Desai, V., von der Linden, A., 

White, S.D.M., 2010, larXiv:1002. 34651 
Pogosyan, D., Pichon, C, Gay, C, Prunet, S., Cardoso, 

J. F., Sousbie, T., & Colombi, S., 2009, MNRAS, 396, 

635 

Porter, S. C, & Raychaudhury, S., 2005, MNRAS, 364, 
1387 

de Putter R., White M., 2005, New Astronomy, 10, 676 
Ragone-Figueroa, Cinthia; Plionis, Manolis, 2007, MN- 
RAS, 377, 1785 
Reblinsky K., Bartelmann M., 1999, A&A, 345, 1 
Sahni, V., Sathyaprakash, B. S., & Shandarin, S. F., 1998, 
ApJL, 495, L5 

Schafer, B. M., 2009, International Journal of Modern 
Physics D, 18, 173 

Schmalzing, J., , 1998, Proceedings of the 12th Pots- 
dam Cosmology Workshop (1997). Eds. V.Mueller, 
S. Gottloeber, J.P.Muecket, J.Wambsganss, 195ff 

Schmalzing, J., Buchert, T., Melott, A. L., Sahni, V., 
Sathyaprakash, B. S., & Shandarin, S. F., 1999, ApJ, 526, 
568 

Shandarin, S.F., Zel'dovich, la. B., 1983, Comments on As- 
trophysics, 10, 33 



Shadarin, S.F., 2004, astro-ph/0405303 



Shandarin, S.F., Sheth, J.V., Sahni, V., 2004, MNRAS, 
353, 162 

Shandarin, S. F., Habib, S., Heitmann, K., 2009, arXiv: 

0912.4471 
Shandarin, S., 2010, larXiv:1011. 19241 
Shaw, L.D., Holder, G.P., Bode, P., 2008, ApJ, 686, 206 
Shen, J., Abel, T., Mo, H. J., & Sheth, R. K., 2006, ApJ, 

645, 783 

Sheth, J. V., Sahni, V., Shandarin, S.F. & Sathyaprakash, 

B., 2003, MNRAS, 343, 22 
De Simone, A., Maggiore, M., Riotto, A., 2010, 

larXiv:1007.1903l 

Skibba, R. A., Sheth, R.K., 2009, MNRAS 392, 1080 
Song, H, Lee, J., 2010, larXiv: 1006.4101 1 
Sousbie, T., Pichon, C, Courtois, H, Colombi, S., & 
Novikov, D., 2008a, ApJL, 672, LI 



© 2010 RAS, MNRAS 000, 000-000 



The geometry of the filamentary environment of galaxy clusters 13 



Sousbie, T., Pichon, C, Colombi, S., Novikov, D., & 

Pogosyan, D., 2008b, MNRAS, 383, 1655 
Sousbie, T., Colombi, S., & Pichon, C, 2009, MNRAS, 393, 

457 

Sousbie, T., 2010. larXiv:1009.4015l 

Sousbie, T., Pichon, C, Kawahara, H., 2010, 
larXiv:1009.4014l 

Stoica, R. S., Martinez, V. J., Mateu, J., & Saar, E., 2005, 

A & A, 434, 423 
Stoica, R. S., M artinez, V. J., & Saar, E., 2008, 

larXiv:0809l358l 

Stoica, R.S., Martinez, V.J., Saar, E., 2009, arXiv: 
0912.2021 

Splinter, R.J., Melott, A. L., Linn, A., Buck, C, Tinker, 

J., 1997, ApJ, 479, 632 
Sunyaev R. A. & Zel'dovich, Ya. B., 1972, 1972, Comm. 

Astrophys, Space Phy, 4, 173 
Sunyaev R. A. & Zel'dovich, Ya. B., 1980, ARA&A, 18, 

537 

Tanaka, M., Finoguenov, A., Kodama, T., Koyama, Y., 

Maughan, B., & Nakata, F., 2009, A & A, 505, L9 
Tormen C, 1997, MNRAS, 290, 411 

van Haarlem, M., van de Weygaert, R., 1993, ApJ, 418, 544 
van Haarlem M.P, Frenk C.S., White S.D.M, 1997, MN- 
RAS, 287, 817 

van de Weygaert, R., Bertschinger, E., 1996, MNRAS, 281, 
84 

van de Weygaert, R., 2002, 2002, Proceedings 2nd Hellenic 

Cos- mology Workshop, 0, 153ff 
van de Weygaert, R., 2006, |arXiv:astro-ph/0607539| 
van de Weygaert, R., & Schaap. W.. 2007. larXiv:0708.1441l 
van de Weygaert, R., Vegter, G., Platen, E., Eldering, B., 

Kruithof, N., 2010. larXiv:1006.2765l 
Voit, CM., 2005, Rev. Mod. Phys. 77, 207 
Wang, H., Mo, H. J., Jing, Y. P., Guo, Y., van den Bosch, 

F. C, & Yang, X., 2009, MNRAS, 394, 398 
Wang, H., Mo, H. J., Jing, Y.P, Yang, X., Wang, Y., 2010, 

larXiv: 1007.0612 1 

Way, M.J., Gazis, PR., Scargle, J.D., 2010, 
larXiv: 1009.03871 

White, R., et al 1999, AJ 118, 2014 
White M., 2001, A&A, 367, 27 

White, M., Hernquist, L., Springel, V., 2002, ApJ, 579, 16 

White M., 2002, ApJS, 143, 241 

White, M., Cohn, J.D., Smit, R., 2010, larXiv: 1005.30221 

MNRAS to appear 
Wojtak, R., Lokas, E.L., Mamon, G.A., Gottlober, S., 

Prada F., Moles M., 2007, A&A 466, 437 
Wray, J.J., Bahcall, N.A., Bode, P., Boettiger, C, Hopkins, 

P.F., 2006, ApJ 652, 907 
Wu, Y., Batuski, D.J., Khalil, A., 2009, ApJ 707, 1160 
Yang X., Mo H.J., van den Bosch F.C., 2008, ApJ, 676, 248 
Zel'dovich, Ia.B., Einasto, J., Shandarin, S.F., 1982, Na- 
ture, 300, 407 

Zhang, Y., Yang, X., Faltenbacher, A., Springel, V., Lin, 
W., Wang, H, 2009, ApJ, 706, 747 



© 2010 RAS, MNRAS 000, 000-000 



